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ABSTRACT 



Aims. The role of carbon monoxide as a cooling agent for the thermal structure of the mid-photospheric to low-chromospheric layers 
of the solar atmosphere in internetwork regions is investigated. 

Methods. The treatment of radiative cooling via spectral lines of carbon monoxide (CO) has been added to the radiation chemo- 
hydrodynamics code C05B0LD. The radiation transport has now been solved in a continuum band with Rosseland mean opacity 
and an additional band with CO opacity. The latter is calculated as a Planck mean over the CO band between 4.4 and 6.2 yum. The 
time-dependent CO number density is derived from the solution of a chemical reaction network. 

Results. The CO opacity indeed causes additional cooling at the fronts of propagating shock waves in the chromosphere. There, the 
time-dependent approach results in a higher CO number density compared to the equilibrium case and hence in a larger net radiative 
cooling rate. The average gas temperature stratification of the model atmosphere, however, is only reduced by roughly 100 K. Also 
the temperature fluctuations and the CO number density are only affected to small extent. A numerical experiment without dynamics 
shows that the CO cooling process works in principle and drives the atmosphere to a cool radiative equilibrium state. At chromospheric 
heights, the radiative relaxation of the atmosphere to a cool state takes several 1000 s. The CO cooling process thus would seem to be 
too slow compared to atmospheric dynamics to be responsible for the very cool temperature regions observed in the solar atmosphere. 
Conclusions. The hydrodynamical timescales in our solar atmosphere model are much too short to allow for the radiative relaxation 
to a cool state, thus suppressing the potential thermal instability due to carbon monoxide as a cooling agent. Apparently, the thermal 
structure and dynamics of the outer model atmosphere are instead determined primarily by shock waves. 

Key words. Sun: chromosphere, photosphere - Hydrodynamics - Radiative transfer - Astrochemistry 



1. Introduction 

• I— I 

, Carbon monoxi de (CO) has been observed in the solar at- 
■ most^ here (e.g.. iNoves & Halll 119721: lAvresI ll98lL lAyres et all 
C3 ' llQSeHSolankiet al.ll 19941; lUitenbroek et al 111994') for more than 
30 years now. The discovery of very low temperatures in ex- 
treme limb observations of strong CO lines at 5 jum by Noyes 
& Hall provoked a crisis in our understanding of the thermal 
structure of the outer layers of the Sun and solar-type stars 
(|Johnsonll973tlWiedemann et alJl994 ). since the very existence 
of CO molecules in significant amounts requires much lower gas 
tem peratures than predic t ed by standard semi -empirical mod- 
els dVernazza etalj|1981l: iFontenla et"anil993l) . Due to its in- 
frared emission, carbon monoxide is known to be a strong cool- 
ing agent that can actively affect the thermal state of the at- 
mospheric gas itself. Therefore, the idea that CO was capable 
of inducing a thermal bifurcation of the solar atmosphere has 
attr acted considerable attention during the past three decades 
(cf.lAvreslll98lHKneeij:1 983; Muchmore & Ulmschneider 1981 
lMuchmorelll986l:ISteffien & Muchmoreil988tlAndersoniil 989). 

Carbon monoxide molecules can form in significant amounts 
under the conditions of the solar photosphere and chromosphere. 
At a given gas density, the equilibrium CO concentration is a 
function of temperature, such that the molecules are quickly 
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dissociated towards higher temperatures. Spectral lines of CO 
at infrared wavelengths can emit radiation effectively and thus 
cool the gas, causing a cooling instability. In the absence of 
compensating mechanical heating, a small perturbation towards 
lower temperature would already allow the formation of more 
CO molecules whose infrared lines would further cool the gas 
and thus induce the formation of evenmoreCO, ultimately lead- 
ing to a "cooling catastrophe" (cf. lAvresll98 l'). This way CO can 
force a cool equilibrium temperature in regions where radiative 
equilibrium (RE) conditions prevail, whereas a critical amount 
of mechanical heating ca n cause a transition to a hot state (e.g., 
[Anderson & Athavll 19891 hereafter AA89). 



In the case of the Sun, the consequences of the 
CO cooling instability have be en investigated by a sim- 
phfied time-dependent approach ( Steffen & Muchmorj 119881 : 
iMuchmore & Uknschneideiill985l: iMuchmore et alj|1988 h that 
treated CO in instantaneous chem ical eq uilibrium (ICE). The 
ICE assumption was also made bv lUitenbroeki (2000) for com- 
puting the spatial distribution of CO in a two-dimensional slice 
from a 3D radiation hydrodynamics simulation. However, this 
assumption is only valid as long as the timescale governing the 
formation of carbon monoxide is short compared to the relevant 
dynamica l time scales. A time-depen dent treatment was pre- 
sented bv lAsensio Ramos et al.l (l2003h . but their simulation was 
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restricted to one spatial dimension in order to keep the problem 
computationally tractable. 

Recently, detailed 2D simulations of the distribution of car- 
bon m onoxide in the solar atmosphere jWedemeyer-Bohm et alj 
I2005L hereafter Paper I) were carried out with the radi- 
ation hydrodyna mics code CO^BOLD (iFreytap et all 120021: 
IWedemeyer et al.l l2004). which can treat chemical reaction net- 
works time-dependently. It was shown that the bulk of CO is 
located in cooler regions in the middle photosphere and that a 
large fraction of all carbon atoms is bound in carbon monoxide 
in the layers above, with the exception of hot propagating shock 
waves (see Fig. [TJ. In contrast to the ICE approach, the high 
shock temperatures do not cause instantaneous destruction of 
CO but rather gradual dissociation of the molecules on chemical 
timescales. Clearly, the highly dynamic and intermittent nature 
of the solar chromosphere makes the time-dependent approach 
mandatory for studying the influence of CO cooling. 

The new two-dimensional simulations presented here were 
done with an upgraded code version that now accounts for 
radiative cooling by CO line s, usin g the simplified treat- 
ment by 'Steff^en & Much mord (1 19881 hereafter SM88) and 
iMuchmo re & Ulmschneiderf ( Il985l liereafter MU85). We de- 
scribe the method and the new simulations in Sects. |2] and |3] 
respectively, and present the results in Sect. |4] followed by dis- 
cussion and conclusions in Sect.|5] 



2. Method 

The applied radiation chemo-hydrodynamics code is described 
in more detail in Paper I. The only difference to the simulations 
presented in Paper 1 concerns the radiative transfer. Instead of 
the grey, i.e., frequency-independent, approach we use two fre- 
quency bands as described by SM88. The first band uses the 
grey Rosseland opacity construct ed from OPA L/PHOENIX 
data dlglesias et all 119921: fHauschild t et al.1 11997|) . It excludes 
the wavelength region around the CO fundamental vibration- 
rotation band in the infra-red beyond 4.6 fim that is accounted 
for in the second band. The opacity in the second band consists 
of the grey Rosseland opacity a:r and an additional Planck mean 
CO opacity kco as a function of gas temperature and CO num- 
ber density. The latter results from the preceding solution of the 
chemical reaction network (see Paper I). 



3. Simulation 

The same initial model as in Paper I is used so that the simu- 
lations with and without CO cooling can be compared directly. 
The model consists of 120 horizontal by 140 vertical grid cells 
with a total extent of 4800 km by 2500 km. It describes a small 
portion of the solar surface layers, just large enough to fit in 
a few granulation cells. The chemical composition is the same 
for each cell in the initial model. The simulation was advanced 
for 86000 s of solar time. The results presented here only refer 
to the last 50000 s, whereas the first 36 000 s are reserved as a 
relaxation phase to ensure the decay of possible initial pertur- 
bations. The lower and upper boundaries are located at heights 
z - -1479 km and z = 1021km, respectively, where the ori- 
gin z - km is defined by the average Rosseland optical depth 
of unity. A correspo nding three-dimensional simulation is cur- 
rently in production dWedemever-Bohm et all in press). 
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Fig. 1. Two-dimensional slices for an exemplary time step at 
t - 70390 s: a) logarithmic gas temperature log T, b) fraction of 
carbon atoms bound in CO nco/«c,totai, c) net radiative heating 
rate Qrad.R for the simulation without CO line cooling, and d) 
the relative difference AQiad between the runs with and without 
CO cooling. Note that the colour range has been limited in the 
last panel in order to to make smaller variations visible. The full 
range is A2,ad - [-2.54,0.20] lO^ergg"' s"'. The lines repre- 
sent contours for T = 5000 K (white solid) and Rosseland opti- 
cal depth unity (long-dashed). The short-dashed line in panels c 
and d marks Qiad.R = and AQiad = 0, respectively. 

4. Results 

4.1. Net radiative lieating rate 

Direct comparison of the two-dimensional simulations with (this 
paper) and without (Paper 1) time-dependent CO opacity shows 
that inclusion of the spectral band of CO causes additional cool- 
ing at the fronts of propagating shock waves in the chromo- 
sphere. There, the time-dependent approach results in a higher 
CO number density compared to the equilibrium case and thus in 
a higher net radiative cooling rate. These small differences make 
the simulations drift apart so that only a statistical comparison 
makes sense. However, we illustrate the effect of CO line cool- 
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Fig. 3. Time evolution towards radiative equilibrium (RE) for 
simulations without (solid) and with CO line cooling for differ- 
ent initial chromospheric temperatures of 2000 K (dot-dashed) 
and 7000 K (triple-dot-dashed): (a) horizontally averaged tem- 
perature stratification for initial time step (thin) and after 
10000 s (thick) and (b) time evolution of the horizontally av- 
eraged temperature at z = 650 km. The shaded area marks the 
transition layer for artificial damping of velocity fluctuations. 



Fig. 2. Horizontally and temporally averaged temperature strati- 
fication (a), temperature fluctuation (b), and CO number density 
(c) for the simulations with (solid line) and without (dashed line) 
radiative CO cooling. The dotted lines represent the difference 
between both cases, e.g. AT" =< Tco >x,t - < Tr >x,i- 



ing by comparing the radiative heating rates calculated with and 
without additional CO line opacity for a given model structure 
(see Fig. [T]). The chosen snapshot is analysed in more detail in 
Paper I. 

For reference, panels a and b of the figure show gas temper- 
ature and the fraction of carbon atoms bound in CO, nco/nc,totab 
respectively. The latter refers to the total number density of all 
carbon atoms, which is here the sum of the densities of included 
species C, CO, and CH. The solid line is a contour for a tem- 
perature of r = 5000 K that outlines the prominent shock wave 
in the low chromosphere and also a high-temperature structure 
belonging to the reversed granulation pattern in the middle pho- 
tosphere. 

The net radiative heating rate 2iad,R for the simulation with- 
out radiative CO cooling is displayed in panel c. The rate repre- 
sents the change in the internal energy of the gas per unit time 
and mass due to radiation. Positive values indicate that energy 
is added, i.e., the gas is heated, whereas negative values indi- 
cate the release of energy via radiative emission, i.e. radiative 
cooling. For most of the computational domain, Qrad.R is close 
to zero or at least small. In contrast, high negative values are 
found in shock waves and close to optical depth unity (the visi- 
ble "surface" of the Sun) where strong radiative emission cools 
the gas. There are also some locations with high positive values, 
indicating radiative heating of the gas. This is obviously the case 
for volume elements that are surrounded by high-temperature re- 
gions, e.g. at X = 3500km, z - -100 km. 

Panel d of Fig.[T]shows the additional contribution of CO to 
the net radiative heating rate, AQrad, which was retrieved by sub- 
tracting the rates from the radiative transfer calculation with and 
without additional CO opacity (Agi-ad = 6rad,co - 6iad,R)- Most 



obvious is the emission (white) at the fronts of the shock waves 
where a significant amount of CO can still eff'ectively cool before 
being dissociated by the approaching hot wave fronts. The effect 
is much smaller outside shock waves. On the other hand, there 
are cool regions with a high CO concentration that absorbs radia- 
tion, resulting in a locally enhanced heating of the gas. Although 
this can make a local difference in the overall contribution of the 
CO lines, which are restricted to a relatively narrow wavelength 
range in the infrared, the overall effect remains minor compared 
to the continuum contribution. 

4.2. Average stratification 

The comparison of the full sequence duration is done by means 
of average stratifications as shown in Fig.|2]for gas temperature, 
root mean square (rms) temperature fluctuation, and CO num- 
ber density. The average temperature of the simulation with CO 
radiative cooling is a bit lower with a maximum deviation of - 
128 K at a height of 495 km, whereas the average temperature 
fluctuation is increased by maximum 124K at z = 105 km. The 
higher fluctuations cause a small increase of 0.25 • 10'^ cm"-' 
(11 %) in average CO number density also at a height of 105 km. 

4.3. Cooling experiment 

The surprisingly small effect of the CO cooling on the thermal 
structure of the atmosphere can be understood from the results 
of a numerical experiment. As the strong dynamics in the atmo- 
sphere disturb the CO formation, we run additional simulations 
in which we damp velocity fluctuations above the middle photo- 
sphere by a "drag-force" characterised by a timescale of Tdamp = 
3 s. After each time step 6t, all spatial velocity components uj 
are reduced with a reduction factor r that ensures a smooth tran- 
sition between the lower undamped domain (z < Z\ - 130 km) 
and the upper damped layers {z> Zi = 310 km): 

vf^a.Jz) = rvf(i,,i,), V; (1) 
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r - I , r > (2) 

"^damp 

^ ^ ziid-zx ^ > ^ > 1 . (3) 

22-^1 

This way the convection in the lower part - 0) remains unaf- 
fected {r - 1), while the chromosphere = 1) can relax to hy- 
drostatic, chemical, and radiative equilibrium. The initial models 
are constructed from the reference model without additional CO 
cooling (see Paper I). The lower part of the model iz < Z\) is 
unchanged. Above we smoothly attach a new atmosphere with 
the gas density set to the horizontal average of the density in the 
original model and the gas temperature changed to a constant 
value of 2000 K or 7000 K (see Fig. |3]a), which are roughly the 
extreme values encountered in the model chromosphere. The ini- 
tial abundances of the molecular species including CO are again 
set to very small numbers. 

Now we run the first simulation with only grey radiative 
transfer with the mean Rosseland opacity band, i.e. frequency- 
independent without CO band, to derive a reference model in 
grey radiative equilibrium. The resulting temperature stratifica- 
tion (see black line in Fig. [3^) reaches an equilibrium value of 
4680 K in the upper photosphere and above. The next simula- 
tions with the additional CO band produce the same tempera- 
ture stratification up to z ~ 150 km. Above that height, CO line 
cooling reduces the mean gas temperature significantly. For in- 
stance, an equilibrium value of ~ 3300 K is reached at a height 
of 500 km, whereas at 650 km the final value is close to 2650 K, 
independent of the initial chromospheric temperature. 

The temporal evolution of the average temperature is shown 
for the chromospheric height of 650 km in panel b. Even for high 
initial temperature, the grey radiative equilibrium is reached af- 
ter only a few minutes. Also the simulation with the cool initial 
chromosphere (2000 K) shows a steeper evolution towards the 
grey RE at the very beginning, but is dominated by the relaxation 
towards the CO-driven cool state afterwards. This second part 
is much slower and approximately exponential with timescales 
of at least ~ 1000 s. The total time for thermal relaxation to a 
cool state is of the order of 3000 s. The hot initial chromosphere 
evolves even slower towards the equilibrium state. In all cases 
the relaxation timescales are long compared to the hydrodynam- 
ical timescales present in the simulations. 

4.4. Prescribed mechanical lieating 

We extended the cooling experiment (see Sect. 14.31 ) and added 
a prescribed (mechanical) heating as done by AA89 (cf. their 
Eq. (1)). The additional heating gadd [erg/g/s] is a function of 
column mass density m : 

Gadd (m) = qo {m/miY e'"''"'" , (4) 

where qo prescribes the heating rate, mo = mi = 0.01 14 and x - 
-0.75. The mean atomic weight yu is set to 2.167 • 10"^^ g/atom. 
The heating term is added to the internal energy s for each grid 
cell for each time step of duration 6t: 

e"™ (/,, 4) = (/,, 4) + 6t gadd iz) . (5) 

The initial model is the same as for the simulation described in 
Sect. [3] The calculations are carried out until the temperature at 
all heights reaches equilibrium. The resulting temperature strati- 
fications are shown in Fig.|4]for different values of qo- For com- 
parison the non-LTE models by AA89 are plotted. We find the 
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Fig. 4. Cooling experiment with prescribed (mechanical) heat- 
ing. Gas temperature after 4h averaged on column mass scale 
for different amounts of prescribed (mechanical) heating with 
log^o - -15.5 and -16.0 (thin solid) and without additional 
heating (thick solid). The thick dot-dashed line represents a sim- 
ulation with grey radiative transfer alone (i.e. without CO). The 
results of AA89 (see their Fig. 1) are plotted for comparison 
(dotted lines). The corresponding mechanical heating (log<7o) is 
noted on the left. T he curve an notated with "RE" is the non-LTE 
RE stratification bv lAndersonI ([T989i) . The shaded area marks the 
transition layer for artificial damping of velocity fluctuations. 

same qualitative dependence of equilibrium temperature stratifi- 
cation on the amount of prescribed heating. Our models have a 
"temperature minimum" in the upper photosphere/low chromo- 
sphere and a steep temperature jump above. Nevertheless there 
are quantitative differences: (i) Our minimum temperatures and 
also the stratifications for no additional heating are cooler than 
those of AA89, (ii) the temperature jump occurs deeper in the 
atmosphere than for AA89, and (iii) the amount of heating re- 
quired for compensating the CO cooling is higher compared to 
AA89. In their calculations the otherwise present cool atmo- 
spheric layer cannot form once the heating exceeds a critical 
value that lies between log^o = -16.0 and -15.5 (see their 
Fig. 1). In contrast, we still get a temperature minimum with 
log^o - -15.5, suggesting that the CO cooling to be compen- 
sated is higher in our models. The different height dependence of 
the average temperature profile is most likely due to differences 
in the treatment of radiative transfer. 

5. Discussion and conclusions 

Our cooling experiment produces a total time for thermal relax- 
ation to a cool state of the order of 3000 s at low-chromospheric 
heights. This result is in line with the earlier results from SM88. 
The final temperature, however, remains ~ 350 K higher in our 
experiment. This can be partly explained by the use of a detailed 
chemical reaction network instead of equilibrium CO densities 
and the corresponding uncertainties in chemical input data in 
both cases. Remaining differences can be attributed to the fact 
that the simulations by Steffen & Muchmore include a station- 
ary flow that differs from the flow field in our simulations and 
the resulting amount of additional cooling due to adiabatic ex- 
pansion. Our cool equilibrium temperature agrees well with the 
model atmospher es in radiative and hydrostatic equilibrium by 
lAndersonI ( Il989h in which CO cooling reduces the tempera- 
ture to 2640 K. It is also in line with AA89, who, depending 
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on the amount of prescribed mechanical heating, find te mpera- 
tures d own to ~ 3000 K, and the 2900 K aheady found bv lAvresI 
(Il981h . 

The treatment of radiative transfer needs to be simpHfied 
in order to keep the problem computationally tractable in the 
framework of a multi-D radiation chemo-hydrodynamics simu- 
lation. The use of only two opacity bands is indeed a strong sim- 
plification, in particular in view of the large number of spectral 
lines of different line strengths in t he CO fundamental vibration- 
rotational bands. l Nordlundl(ll985h states that "the cooling in the 
infrared CO lines is compensated for by a heating in the ODFs 
in the blue and near UV" if UV opacities (A < 400 nm) are 
treated as pure absorption, as in the method applied here (based 
on MU85). T he alternative as sumption of pure scattering is ques- 
tionable (cf. lAndersonlll989i) . The fact that we still find a net 
cooling due to CO in the absence of atmospheric dynamics can 
be interpreted such that our simplified LTE tr eatment overesti- 
mates the cooling ability of the IR CO lines. iMuchmore et alj 
(119881) compare the scheme by MU85 with a more detailed 
scheme with more frequency points and conclude that (i) the CO 
overtone bands at 2.2 /im play only a minor role compared to 
the fundamental bands in the 5 /vm region and (ii) the CO cool- 
ing might be overestimated by a factor of three with this simple 
scheme. 

A similar conclusion can be drawn from the comparison of 
the temperature stratifications from our extended cooling exper- 
iment with the NLTE models by AA89 (see Sect. l4!4l l. It re- 
veals that the equilibrium temperature in the atmosphere is lower 
than found by AA89 and that the amount of additional heating 
required to compensate CO cooling is higher for our models. 
Obviously, our approach represents an upper limit for the net 
radiative cooling rate. So if even an over-efficient CO cooling 
mechanism cannot induce significant effects on the atmospheric 
stratification, as shown in the study presented here, then we can 
safely conclude that a more realistic treatment with smaller net 
cooling will have even less effect. 

The inclusion of CO line opacity in our simulations in- 
deed produced additional radiative cooling at the front of fast- 
propagating shock waves in the upper atmospheric layers due 
to a higher CO number density compared to the instantaneous 
chemical equilibrium (ICE) case. In contrast to the ICE case, 
the dissociation of CO proceeds on a finite timescale so that the 
CO number density in the front of a passing shock wave is re- 
duced gradually and not instantaneously. As a consequence, ad- 
ditional radiative cooling is present at these locations. The re- 
sulting changes in the thermal structure, however, remain small. 

It is commonly known that the dynamics tend to be too 
strong in 2D compared to 3D. But even in our 3D simulations 
we do not find the passage of shock waves to be less frequent, 
and the temperature fluctuations are comparable in amplitude ex- 
cept for somewhat larger differences i n the middle photosphere 
(Wedemever 2003). The 3D model by IWedemever-Bohm et 
din pressi) nevertheless yields a CO distribution that is very simi- 
lar to the 2D results presented here, confirming that the dynamics 
in the 2D model are reasonably realistic. 

We conclude that the self-amplifying cooling process of CO 
can in principle operate in the solar atmosphere (as demonstrated 
in Fig. [3] and by earlier works) and could lead to a bifurcation 
between cool areas dominated by RE and warmer areas dom- 
inated by mechanical heating. However, a relaxation to a cool 
state of the solar atmosphere is prevented in the presence of 
pronounced atmospheric dynamics as predicted by our models. 
The crucial point is that the CO radiative relaxation timescales 
are too long compared to the frequent passage of shock waves. 



Instead, the co-existence of hot and cool regions is predomi- 
nantly caused by mechanical heating due to these propagating 
shock waves and adiabatic cool ing of the resulting post-shock re- 
gions jWedemever et al.ll2004l) . Without any doubt, the present 
state-of-the-art simulations still admit to some limitations, and 
the predictions of the velocity and temperature fluctuations in 
the higher layers are thus affected by uncertainties. Detailed ob- 
servational tests of atmospheric properties, such as temporal and 
spatial temperature and velocity variations and their centre-to- 
limb behaviour, are therefore highly desired. 
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